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Abstract 



In this paper, we complete the global qualitative analysis of the well-known Fitz- 
Hugh-Nagumo neuronal model. In particular, studying global limit cycle bifurca- 
. tions and applying the Wintner-Perko termination principle for multiple limit cy- 

• I cles, we prove that the corresponding dynamical system has at most two limit cycles. 
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. 1 Introduction 

o ^ 

We consider tlie well-known FitzHugli-Nagumo model in the form 



W = e{V-6W), 

where V is the membrane potential, Ty is a recovery variable, and / is the 
magnitude of stimulus current, which is a two-dimensional simplification of 
the classical Hodgkin-Huxley model of the spike dynamics in a biological neu- 
ron [5], [H], [IS], [H], [20]. This system was suggested by FitzHugh (1961) [5], 
who called it "Bonhoeffer - van der Pol model", and the equivalent circuit was 
constructed by Nagumo et al. (1962) [16]. The motivation for the FitzHugh- 
Nagumo model was to isolate conceptually the essentially mathematical pro- 
perties of excitation and propagation from the electrochemical properties of 
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sodium and potassium ion flow. The model consists of a voltage-like variable 
having cubic nonlinearity that allows regenerative self-excitation via a posi- 
tive feedback, and a recovery variable having a linear dynamics that provides 
a slower negative feedback. While the Hodgkin-Huxley model is more real- 
istic and biophysically sound, only projections of its four-dimensional phase 
trajectories can be observed. The simplicity of the FitzHugh-Nagumo model 
permits the entire solution to be viewed at once. This allows a geometrical ex- 
planation of important biological phenomena related to neuronal excitability 
and spike-generating mechanism [T3]. 

The phase portrait of the FitzHugh-Nagumo model (1.1) depicts the ^-nuU- 
cline, which is the A^-shaped curve obtained from the condition V = 0, and 
the ly-nullcline, which is a straight line obtained from the condition W = 0. 
The intersection of nuUclines is an equilibrium (a singular point) of the system 
(1.1), which may be unstable if it is on the middle branch of the \^-nullcline, 
i. e., when I is strong enough. In this case, the model exhibits periodic (tonic 
spiking) activity. 

The FitzHugh-Nagumo model explained the absence of all-or-none spikes in 
the Hodgkin-Huxley model in response to stimuli, i.e., pulses of the injected 
current I. Weak stimuli (small pulses of I) result in small- amplitude tra- 
jectories that correspond to subthreshold responses; stronger stimuli result 
in intermediate-amplitude trajectories that correspond to partial-amplitude 
spikes; and strong stimuli result in large-amplitude trajectories that corre- 
spond to suprathreshold response — firing a spike. 

Similarly to the Hodgkin-Huxley model, the FitzHugh-Nagumo model does 
not have a well-defined firing threshold in the absence of a saddle equilibrium. 
This feature is the consequence of the absence of all-or-none responses. The ap- 
parent illusion of threshold dynamics and all-or-none responses in both models 
is due to the existence of the "quasi-threshold" , which is a canard trajectory 
that follows the unstable (middle) branch of the A^-shaped l^-nullcline. 

The FitzHugh-Nagumo model explains the excitation block phenomenon, i. e., 
the cessation of repetitive spiking as the amplitude of the stimulus current in- 
creases. When / is weak or zero, the equilibrium (intersection of nullclines) is 
on the left (stable) branch of l^-nuUcline, and the model is resting. Increas- 
ing / shifts the nuUcline upward and the equilibrium slides onto the middle 
(unstable) branch of the nuUcline. The model exhibits periodic spiking ac- 
tivity in this case. Increasing the stimulus further shifts the equilibrium to 
the right (stable) branch of the iV-shaped nuUcline, and the oscillations are 
blocked (by excitation). The precise mathematical mechanism involves ap- 
pearance and disappearance of a limit cycle attractor, and it is reviewed in 
detaU by Izhikevich (2007) [H]. 
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This model explained also the phenomenon of post-inhibitory (rebound) spikes, 
called anodal break excitation at that time. As the stimulus I becomes neg- 
ative (hyperpolarization), the resting state shifts to the left. As the system 
is released from hyperpolarization (anodal break), the trajectory starts from 
a point far below the resting state (outside the quasi-threshold) , makes a 
large-amplitude excursion, i. e., fires a transient spike, and then returns to the 
resting state. 

The FitzHugh-Nagumo model explained the dynamical mechanism of spike 
accommodation in Hodgkin-Huxley-type models. When stimulation strength 
increases slowly, the neuron remains quiescent. The resting equilibrium of (1.1) 
shifts slowly to the right, and the state of the system follows it smoothly with- 
out firing spikes. In contrast, when the stimulation is increased abruptly, even 
by a smaller amount, the trajectory could not go directly to the new resting 
state, but fires a transient spike; see figure. Geometrically, this phenomenon 
is similar to the post-inhibitory (rebound) response. 

The FitzHugh-Nagumo equations became a favorite model for reaction-diffusion 

systems 

V = I _w -aV + ia + 1)^2 _ y3 ^ Y 

^ ' (1.2) 

W = e{V - 6W), 

which simulate propagation of waves in excitable media, such as heart tissue or 
nerve fiber. Here, the diffusion term Vxx is the second derivative with respect 
to the spatial variable x. Its success is mostly due to the fact that the model 
is analytically tractable, and hence it allows derivation of many important 
properties of traveling pulses without resort to computer simulations. 

Without loss of generality, the system (1.1) can be written in the canonical 
form 

X = (7 (5 — 1) ?/ + (7 — a) X -|- fex^ — cx^ = P(x, y), 

y = x-5y = Q{x,y). 

Such a system was studied earlier, e. g., in [T9]. However, its qualitative analy- 
sis was incomplete, since the global bifurcations of multiple limit cycles could 
not be studied properly by means of the methods and techniques which were 
used earlier in the qualitative theory of dynamical systems. Applying new 
bifurcation methods and geometric approaches developed in [1], [6]-[13], we 
complete the qualitative analysis of the FitzHugh-Nagumo model and prove, 
in particular, that the corresponding dynamical system (1.3) has at most two 
limit cycles. In Sections 2-4 of this paper, we recall basic facts and results 
from the global bifurcation theory of polynomial dynamical systems and its 
applications. These results, together with the methods of [3], [1], [6]-[T3], are 
used in Sections 5, 6 for the study of singular point and limit cycle bifurcations 
of the system (1.3). 
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2 Preliminaries 



In this paper, geometric aspects of Bifurcation and Catastrophe Theories are 
used and developed [6], |T5], [18]. First of all, the two-isocline method which 
was developed by Erugin is used, see [6]. An isocline portrait is the most na- 
tural construction for a polynomial equation. It is sufficient to have only two 
nullclines (or isoclines of zero and infinity in our terminology) to obtain princi- 
pal information on the original polynomial system, because these two isoclines 
are right-hand sides of the system. Geometric properties of isoclines (conies, 
cubics, quartics, etc.) are well-known, and all isocline portraits can be easily 
constructed. By means of them, all topologically different qualitative pictures 
of integral curves to within a number of limit cycles and distinguishing center 
and focus can be obtained. Thus, it is possible to carry out a rough topological 
classification of the phase portraits for the polynomial dynamical systems. It 
is the first application of Erugin's method. After studying contact and rota- 
tion properties of the isoclines, the simplest (canonical) systems containing 
limit cycles can be also constructed. Two groups of parameters can be dis- 
tinguished in such systems: static and dynamic. Static parameters determine 
the behavior of phase trajectories in principle, since they control the number, 
position, and character of singular points in a finite part of the plane (finite 
singularities). The parameters from the first group determine also a possible 
behavior of separatrices and singular points at infinity (infinite singularities) 
under variation of the parameters from the second group. The dynamic pa- 
rameters are field rotation parameters, see [T], [B], [IB]- They do not change 
the number, position and index of the finite singularities, but only involve the 
vector field in a directional rotation. The rotation parameters allow to control 
the infinite singularities, the behavior of limit cycles and separatrices. The 
cyclicity of singular points and separatrix cycles, the behavior of semi-stable 
and other multiple limit cycles are controlled by these parameters as well. 
Therefore, by means of the rotation parameters, it is possible to control all 
limit cycle bifurcations and to solve the most complicated problems of the 
qualitative theory of dynamical systems. 

In [S], [7], [H], [in], [13], some complete results on quadratic systems have 
been presented. In particular, it has been proved that for quadratic systems 
four is really the maximum number of limit cycles and (3 : 1), i.e., three 
limit cycles around one focus and the only limit cycle around another focus, 
is their only possible distribution (this is a solution of Hilbert's Sixteenth 
Problem in the quadratic case of polynomial dynamical systems). In [S], some 
preliminary results on generalizing new ideas and methods of [0] to cubic 
dynamical systems have already been established. In particular, a canonical 
cubic system of Kukles type has been constructed and the global qualitative 
analysis of its special case corresponding to a generalized Lienard equation has 
been carried out. It has been proved also that the foci of such a Lienard system 
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can be at most of second order and that such system can have at most three 
hmit cycles on the whole phase plane. Moreover, unlike all previous works 
on the Kukles-type systems, global bifurcations of limit and separatrix cycles 
using arbitrary (including as large as possible) field rotation parameters of the 
canonical system have been studied in [8|. As a result, the classification of all 
possible types of separatrix cycles for the generahzed Lienard system has been 
obtained and all possible distributions of its limit cycles have been found. In 
|12j . a solution of Smale's Thirteenth Problem proving that the Lienard 
system with a polynomial of degree 2k + 1 can have at most k limit cycles 
has been presented. In [4^, we have completed the global qualitative analysis 
of a quartic ecological model. All of these methods and results can be applied 
to the global qualitative analysis of the FitzHugh-Nagumo neuronal model 
as well. 

In [3], we have already carried out the global qualitative analysis of a poly- 
nomial dynamical system as a learning model of neural networks [2], [T7] . 
Learning models are algorithms, implementable as neural networks, that aim 
to mimic an adaptive procedure. A neural network is a device consisting on 
interconnected processing units, designated neurons. An input presented to 
the network is translated as a numerical assignment to each neuron. This will 
create a sequence of internal adjustments leading to a learning process. An 
input vector, denoted by C,, represents an n-dimensional random vector with 
independent components. This means that the joint probability distribution 
function is the product of n density functions. The output value, denoted by V, 
is the outcome of the network's action on ^ and is given by Y^^=i^j^j-, where 
ujj is the connecting weight for the synapse attached to the input neuron j. 
Since new synapses may be created under a constant error rate, E, a synaptic 
strength may capture nearby activity. This is done by the creation of tempo- 
rary synapses from the closest neurons to the output one. The synaptic rate 
of change is given by 

= V{{1 - E)i, + (i?/2)(e.+i + ^.-i) - Vui), 

for 2 = 2, . . . , n — 1, or 

uj, = V{{1 - E)i, + (E/2)e.±i - Vuji), 

for i = 1 or n, respectively. Substituting the value of V in the expression of 
Cji, we obtain 

z 7^ 1 and n, 



= < 



z = 1 or n, 
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what can be reduced to the equation 

^ = TClo-(lo,Cuj)uj, (2.1) 
at 

with C = = ^ = {^1,^2, • • • , ^n}, a correlation matrix of ex- 

pected values, and T, a tridiagonal substochastic matrix given by tij = if 
> l,tij = E/2 ii = 1, and ta = 1 — E, for all i and j = 1, . . . , n 

Desirable initial conditions are those with trajectories that converge to some 
equilibrium (singular) point of (2.1). This will assume a natural weight assign- 
ment as a result of the learning process. Knowledge on the stability of equilibria 
provides information on the robustness of the learning process. Existence of 
cycles might represent a different kind of learning where a whole continuum 
of connecting weight vectors emerges instead of just a single vector. In [3J, we 
have restricted our attention to two dimensions. 

For two input neurons, (2.1) can be written as a cubic dynamical system 
X = {{l-e)a + {e/2)b)x + {{l-e)b+{e/2)c)y - x{ax'^ + 2bxy+cy'^), 

(2.2) 

y = {{e/2)a+{l-e)b)x + {{e/2)b+{l-e)c)y - y{ax'^ + 2bxy + cy'^), 

where the parameters e and a, b, c represent, respectively, the probability of 
synaptic formation and the weight strengths for the synapses attached to the 
input neurons [3]. Thus, we have got a four-parameter planar dynamical sys- 
tem for investigation. Applying techniques based both on classical Poincare 
and Dulac methods and also on some methods developed in [6] , we have stud- 
ied the global bifurcations of singular points and limit cycles of the cubic 
system (2.2), a learning model of planar neural networks [3]. 

Some of these techniques can be extended to higher- dimensional dynamical 
systems [6], [15], [18]. So, for the global analysis of limit cycle bifurcations 
(in particular, for solving the uniqueness problem) we have used the Perko 
planar termination principle stating that the maximal one-parameter family 
of multiple limit cycles terminates either at a singular point, which is typi- 
cally of the same multiplicity, or on a separatrix cycle, which is also typically 
of the same multiplicity [18j. This principle is a consequence of the Wintner 
principle of natural termination, which was stated for higher-dimensional dy- 
namical systems (see [H]), where one-parameter families of periodic orbits 
of the restricted three-body problem are studied and Puiseux series are used 
to show that in the analytic case any one-parameter family of periodic orbits 
can be uniquely continued through any bifurcation except a period-doubling 
bifurcation. Thus, the Wintner-Perko termination principle and the method 
developed in [3], [1], [6]-[T3] can be applied to the further global qualitative 
analysis of neural dynamical systems. 
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3 Basic facts on limit cycles 



Consider a polynomial dynamical system in the vector form 

x = f{x,n), (3.1) 

where xeR^, ^ G R"; / G ( / is a polynomial vector function). 

Let us recall some basic facts concerning limit cycles of (3.1). But first let us 
state two fundamental theorems from the theory of analytic functions [B]. 

Theorem 3.1 (Weierstrass Preparation Theorem). Let F{w,z) be an 
analytic in the neighborhood of the point (0, 0) function satisfying the following 
conditions 

aw a w a w 



Then in some neighborhood \w\ < e, \z\ < 6 of the points (0,0) the function 
F[w, z) can be represented as 

F{w, z) = {w'' + Ai{z)w''~^ + ... + Ak-iiz)w + Ak{z))^{w, z), 

where ^{w, z) is an analytic function not equal to zero in the chosen neigh- 
borhood and Ai{z), . . . , Ajs{z) are analytic functions for \z\ < 6. 

From this theorem it follows that the equation F{w, 2;) = in a sufficiently 
small neighborhood of the point (0, 0) is equivalent to the equation 

+ Ai{z)w''^^ + ... + Ak^i{z)w + Ak{z) = 0, 

which left-hand side is a polynomial with respect to w. Thus, the Weierstrass 
preparation theorem reduces the local study of the general case of implicit 
function w{z), defined by the equation F{w,z) = 0, to the case of implicit 
function, defined by the algebraic equation with respect to w. 

Theorem 3.2 (Implicit Function Theorem). Let F[w,z) be an analytic 
function in the neighborhood of the point (0, 0) and F{0, 0) =0, -F^(0, 0) 7^0. 

Then there exist 6 > and e > such that for any z satisfying the condition 
\z\ < 6 the equation F{w, z) = has the only solution w = f{z) satisfying the 
condition \ f{z)\ < e. The function f{z) is expanded into the series on positive 
integer powers of z which converges for \z\ < 6, i. e., it is a single-valued 
analytic function of z which vanishes at z = 0. 
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Assume that the system (3.1) has a hmit cycle 

of minimal period Tq at some parameter value fi = fiQ E R" (Fig. 1). 




1 2 



FIG. 1. The Poincare return map in the neighborhood of a muhiple hmit cycle. 

Let I be the straight line normal to Lq at the point Pq = <Pq{0) and s be 
the coordinate along I with s positive exterior of Lq. It then follows from the 
implicit function theorem that there is a 5 > such that the Poincare map 
h{s,i-i) is defined and analytic for |s| < 6 and \\fj, — /Xq|| < 6. Besides, the 
displacement function for the system (3.1) along the normal line / to Lq is 
defined as the function 

d{s, n) = h{s, /x) — s. 

In terms of the displacement function, a multiple limit cycle can be defined as 
follows |6]. 

Definition 3.1. A limit cycle Lq of (3.1) is a multiple limit cycle iff d{0, /Xq) = 
dr{0,i-iQ) = and it is a simple limit cycle (or hyperbolic limit cycle) if it is 
not a multiple limit cycle; furthermore, Lq is a limit cycle of multiplicity m iff 

rf(0,^o) = ^.(0, A^o) = • • • = d^-'\0,fio) = 0, 4'"Ho,/Xo) ^ 0. 

Note that the multiplicity of Lq is independent of the point Pg ^ -^o through 
which we take the normal line /. 

Let us write down also the following formulas which have already become 
classical ones and determine the derivatives of the displacement function in 
terms of integrals of the vector field / along the periodic orbit (foit) |6]: 

To 

rf,(0,/Xo) = exp y'v-/((^oW>^o)dt-l (3.2) 



and 
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To / t 



Jexp l-Jv- fi^oir), Mo) dr / A f^^i^oit), Mo) (3.3) 



ll/(¥'o(0),/Xo)|| ^ , ^ 

for j = l,...,n, where cjo = ±1 according to whether Lq is positively or 
negatively oriented, respectively, and where the wedge product of two vectors 
X = (xi, X2) and y = {yi, 2/2) in R-^ is defined as 

X Ay = xiy2 - X2 yi- 



Similar formulas for dss{0, fi^) and ds^-{0, /-Iq) can be derived in terms of in- 
tegrals of the vector field / and its first and second partial derivatives along 
^o{t). The hypotheses of theorems in the next section will be stated in terms 
of conditions on the displacement function d{s, fi) and its partial derivatives 
at (0,Mo) [6]. 



4 Bifurcation surfaces of multiple limit cycles 



In this section, we restate Perko's theorems on the local existence of (n— m+1)- 
dimensional surfaces. Cm, of multiplicity-m limit cycles for the polynomial sys- 
tem (3.1) with fi G R" and n > m > 2. These results describe the topological 
structure of the codimension (m — 1) bifurcation surfaces Cm- For m = 2,3,4, 
C2, C3, and C4 are the familiar fold, cusp, and swallow-tail bifurcation sur- 
faces; for m > 5, the topological structure of the surfaces Cm is more complex. 
For instance, C5 and Cq are the butterfly and wigwam bifurcation surfaces, 
respectively [T8]. Since the proofs of the theorems in this section, describing 
the universal unfolding near a multiple limit cycles of (3.1), parallel the clas- 
sical proofs of Catastrophe Theory, we will only state the theorems (see [IB] 
for more detail). 

Definition 4.1. An (n— l)-dimensional analytic surface C2 C R" is an [n—l)- 
dimensional fold bifurcation surface of multiplicity-two limit cycles of (3.1) 
through a point /Xg G -R", if for all e > there exists a 5 > such that for each 
G C2 with 11^ — ^igll < (5, the system (3.1) has a unique multiplicity-two 
limit cycle in an e-neighborhood of Lq and the system (3.1) undergoes a 
fold bifurcation at L^; i.e., for \\^ — /Xq|| < 5, splits into a simple stable 
and a simple unstable limit cycles in an e-neighborhood of Lq for /x on one 
side of C2 and L^ vanishes for ^ on the other side of C2. Cf. Fig. 2. 

Theorem 4.1. Suppose thatn > 2, that for fi = /Xq G -R" the system (3.1) has 
a multiplicity-two limit cycle Lq, and that d^-^{0, fi^) 7^ 0. Then given e > 0, 

there is a 6 > and a unique function g{n2-, • • • , A^n) with 5'(/i2°\ • • • > A^n^'') = 
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FIG. 2. The fold bifurcation surface. 

fif\ defined and analytic for — f^'i^l < ■ ■ ■ , |/^n — < such that for 
\l^2 - 122^1 < 6,..., - ) I < 5, 

C2 : /ii = 5'(/i2, . . . 

is an [n — 1)- dimensional, analytic fold bifurcation surface of multiplicity-two 
limit cycles of (3.1) through the point ^q. 

Definition 4.2. An analytic surface C3 C R" is an {n — 2) -dimensional cusp 
bifurcation surface of multiplicity-three limit cycles of (3.1) through a point 
Hq G i?"", if for all e > there exists a 6 > such that for each ^ E with 
11/^ — /XqII < 5, the system (3.1) has a unique multiplicity-three limit cycle 
in an e- neighborhood of Lq and the system (3.1) undergoes a cusp bifurcation 
at L^; i.e., C3 is the intersection of two (n— l)-dimensional fold bifurcation 
surfaces of multiplicity-two limit cycles of (3.1), C^, which intersect in a cusp 
along C3; for ||/^ — /Xq|| < 5 and for ^ in the cuspidal region between C2 
and C2 (shaded in Fig. 3), the system (3.1) has three simple limit cycles in 
an e- neighborhood of Lq; and for — /Xq|| < 5 and ^ outside the cuspidal 
region, the system (3.1) has one simple limit cycle in an e-neighborhood of Lq. 
Cf. Fig. 3. 

Theorem 4.2. Suppose that n > 3, that for n = E i2" the system (3.1) 
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FIG. 3. The cusp bifurcation surface. 



has a multiplicity-three limit cycle Lq, that (i^i(0, /Xq) 0; drij.i{0, Mo) '^'^^ 
forj = 2, . . . ,n, 

d{d, dr 



Then given e > 0, there is a 6 > and constants Uj = ±1 forj = 2, . . . , n, and 
there exist unique functions hi{fi2, . . . , fin), ^2(/^25 ■ ■ ■ , yUn) (iiT'd g^{fi2, ■ ■ ■ , fJ^n) 
with h,{f,^^\ . . . = /if), . . . ) = /if) and ^/±(/if\ . . . ,/iW) = 

/if), where hi and h2 are defined and analytic for |/ij — /if ) | < 5, j = 2, . . . , ra. 



and are defined and continuous for < crj(/ij — /if )) < 5 and analytic f 



or 



< ^^(/Xj — /if)) < 5, j = 2, . . . , n siic/i t/iat 



: 



/il = /il(/i2, . . .,IJ,n) 
/il = /i2(/i2, • • • ,/in) 



is an (n, — 2) -dimensional, analytic, cusp bifurcation surface of multiplicity- 
three limit cycles of (3.1) through the point and 

Ct : /il = 5'=^(/i2, . . . ,/i„) 

are two {n — 1)- dimensional, analytic, fold bifurcation surfaces of multiplicity- 
two limit cycles of (3.1) which intersect in a cusp along C3. 
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Definition 4.3. An analytic surface C4 C R'^ is an {n — 3) -dimensional 
swallow-tail bifurcation surface of multiplicity-four limit cycles of (3.1) through 
a point /Xq G -R", if for all e > there exists a 5 > such that for each /x G C4 
with ll/x — /XqII < 6, the system (3.1) has a unique multiplicity- four limit cycle 
in an e- neighborhood of Lq and the system (3.1) undergoes a swallow-tail 
bifurcation at L^; i.e., C4 is the intersection of two (n — 2)-dimensional cusp 
bifurcation surfaces of multiplicity-three limit cycles which intersect in a 
cusp along C4; furthermore, there are three (ra— l)-dimensional fold bifurcation 
surfaces of multiphcity-two limit cycles of (3.1), C2*'', « = 0, 1, 2, such that C2°'* 
and €2^^ intersect in a cusp along C^, C2°'* and 6*2^-' intersect in a cusp along 
, and Cq^^ and €2^^ intersect along an (n— 2)-dimensional surface on which 
(3.1) has two multiplicity-two limit cycles; finally, for — /XqH < S and for /x 
in the swallow-tail region (shaded in Fig. 4), the system (3.1) has four simple 
limit cycles in an e-neighborhood of Lq; for — /Xq|| < S and fx above the 
surfaces 6*2*^ i = 0,1,2, the system (3.1) has two simple limit cycles in an 
e-neighborhood of Lq; and for ||/x — ^tQ|| < 6 and /x below the surfaces C2*'*, 
i = 0,1,2, the system (3.1) has no limit cycles in an e-neighborhood of Lq. 
Cf. Fig. 4. 




FIG. 4. The swallow-tail bifurcation surface. 



Theorem 4.3. Suppose that n > 4, that for ^ = /Xq G -R" the system (3.1) 
has a multiplicity-four limit cycle Lq, that (i^i(0,^o) 7^ 0, dr^^{0, Hq) 7^ 0, 
drrfii (0, /Xq) 7^ 0, and that for j = 2, . . . ,n, 

®'''*'-(0.M„)/0. |^(0.M„)^0. t?^(O.M„)^0. 



a(/ii,/ij) ' ' <9(/ii,/Xj) ' ' 9(/ii,/ij) 



12 



Then given e > 0, there is a S > and constants Ujk = ±1 for j = 2, . . . ,n, 
= 1, 2, and there exist unique functions 5'i(/i2, • • • , /in), • • • , /^n) o-nd 

Fi(/i2, ...,Hn), with gi{n^2\ /^n ^) = h^{f^2\ ■■■,f^n^) = Fi{nf\ fi^^'^) = 

IJ,i'\ for i = 0,1,2 and k = 1,2, where Fi is defined and analytic for i = 
0, 1, 2, and \jjj — jjf^l < 6, j = 2, . . . ,n, are defined and continuous for 

< oJjkifJ'j — fif^) < 6 and analytic for < u)jk{fJ'j — f^f^) < ^, j = 2, . . . , n, 
k = 1,2, and for i = 0, 1,2, Qi is defined and analytic in the cuspidal region 
between the surfaces /xi = hf{n2, ■ ■ ■ which intersect in a cusp, and gt is 
continuous in the closure of that region, such that 









■ ■ 5 f^n) 


C4 : < 


/ii 


= i^i(/i2,. 






/ii 







is an {n — 3) -dimensional, analytic, swallow-tail bifurcation surface of multi- 
plicity-four limit cycles of (3.1) through the point /-i^ which is the intersection 
of two {n — 2) -dimensional, analytic, cusp bifurcation surfaces of multiplicity- 
three limit cycles of (3.1 ), 

A^i = hf{^i2, • • • 

/il = /if (/i2, ■ • ■ 

which intersect in a cusp along C4; furthermore, = c!f^ f] ^^2^^ f*^'^ = 
Cf'^ n Cf^ where for i = f), 1,2, 

fil = gi{fi2, ■ ■ ■ , fJ^n) 

are {n — 1)- dimensional, analytic, fold bifurcation surfaces of multiplicity-two 
limit cycles of (3.1) which intersect in cusps along and in an [n — 2)- 
dimensional, analytic surface C'2'' f] Cf"^ on which (3.1) has two multiplicity- 
two limit cycles (Fig. 4 and Fig. 5). 



Based on Theorems 3.1, 3.2, the following generalization of Theorems 4.1-4.3 
can be proved on induction [T8] . 



Theorem 4.4. Given m >2. Suppose that n > m, that for = /Xq G -R"" the 
polynomial system (3.1) has a multiplicity-m limit cycle Lq, that 

^^(0,Mo)7^0, 1^(0,^0)7^0,..., (0,/Xo)7^0, 



ojjLi OjjLi a/ii 
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Q (Q) 




FIG. 5. The bifurcation curve (one-parameter family) of multiple limit cycles. 



for i, j = 0, . . . ,m — 2 with i ^ j and k = 2, . . . ,n. 

Then given e > there is a S > such that for W/j, — HqW < S, the system (3.1) 
has 

(1) a unique {n — m + 1) -dimensional analytic surface Cm of multiplicity-m 
limit cycles of (3.1) through the point Hq] 

(2) two {n — m + 2) -dimensional analytic surfaces Cm-i of multiplicity-{m—l) 
limit cycles of (3.1) through the point ^tg which intersect in a cusp along Cm', 

(j) exactly j, {n—m-\-j)-dimensional analytic surfaces Cm-j+i of multiplicity- 
(m — j + 1) limit cycles of (3.1) through the point which intersect pairwise 
in cusps along the bifurcation surfaces Cm-j+2', 

(m — 1) exactly (m — 1), {n — 1) -dimensional analytic fold bifurcation surfaces 
C2 of multiplicity-two limit cycles of (3.1) through the point (Xq which intersect 
pairwise in a cusp along the {n — 2) -dimensional cusp bifurcation surfaces C3. 



and that 
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5 Singular points of the FitzHugh— Nagumo system 



The study of singular point of the system (1.3) will use two index theorems 
by H. Poincare, see |T]. But first let us define the Poincare index [1]. 

Definition 5.1. Let S* be a simple closed curve in the phase plane not passing 
through a singular point of the system 



where P{x, y) and Q{x, y) are continuous functions (for example, polynomials), 
and M be some point on S. If the point M goes around the curve S in positive 
direction (counterclockwise) one time, then the vector coinciding with the 
direction of a tangent to the trajectory passing through the point M is rotated 
through the angle 27rj (j = 0, ±1, ±2, . . .). The integer j is called the Poincare 
index of the closed curve S relative to the vector field of system (5.1) and has 
the expression 



According to this definition, the index of a node or a focus, or a center is equal 
to +1 and the index of a saddle is —1. 

Theorem 5.1 (First Poincare Index Theorem). If Nf, Nc, and C are 

respectively the number of nodes, foci, centers, and saddles in a finite part of 
the phase plane and N' and C are the number of nodes and saddles at infinity, 
then it is valid the formula 



Theorem 5.2 (Second Poincare Index Theorem). If all singular points 
are simple, then along an isocline without multiple points lying in a Poincare 
hemisphere which is obtained by a stereographic projection of the phase plane, 
the singular points are distributed so that a saddle is followed by a node or a 
focus, or a center and vice versa. If two points are separated by the equator 
of the Poincare sphere, then a saddle will be followed by a saddle again and a 
node or a focus, or a center will be followed by a node or a focus, or a center. 

Consider the system (1.3). Its finite singularities are determined by the alge- 
braic system 



x = P{x,y), y = Q{x,y) 



(5.1) 




N + Nf + N, + N' = C + C' + l. 



{'y 5 — 1) y + — a) X + bx"^ 
X — 6 y = . 



— cx^ = 0, 



(5.2) 
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From (5.2) [T], [6], [19], we will get a singular point (0,0) of antisaddle-type 
(i.e., a node, a focus, or a center) and at most two points (a saddle and an 
antisaddle or, if one, a saddle-node) defined by the condition 

2 , 7 (5 — 1 X 
cx —ox 7 + a = 0, y = ■ (5.3) 



To get singular points at infinity, consider the corresponding differential equa- 
tion 

dy X — 6y 



dx {'J 6 — 1) y + {'-y — a) X + bx'^ — cx^ 



(5.4) 



Dividing the numerator and denominator of the right-hand side of (5.4) first 
by x^ {x 7^ 0), denoting y/x by u, and then by y^ {y ^ 0), denoting xjy by w, 
we will get two infinite singularities: u = (a simple node in the direction of 
the X-axis) and f ^ = (a triple saddle in the direction of the y-axis), see PP, [S]- 



6 Global limit cycle bifurcations in the system 



To investigate global limit cycle bifurcations in the system (1.3), we will use 
the results of the previous sections and will apply the method developed in [3], 
[1], [6]-[l3]. The sense of this method is to obtain the simplest (well-known) 
system by vanishing some parameters (usually field rotation parameters) of 
the original system and then to input these parameters successively one by 
one studying the dynamics of limit cycles on the whole phase plane. 

Let us study rotation properties of the parameters of (1.3). Applying the 
definition of a field rotation parameter (i.e., a parameter which rotates the 
field in one direction [1], |6], [H]), to the system (1.3) written in the form 

a; = i?(x,y) + 7(5(x,?/) = P(x,y), y = Q{x,y), (6.1) 

where R{x,y) = —y — ax + bx"^ — cx^ and Q{x,y) = x — 6y, calculate the 
corresponding determinant for the parameter 7: 

= pg; - gp; = -g^ < o. (6.2) 

It follows from (6.2) that on increasing 7 the vector field of (1.3) is rotated in 
negative direction (clockwise) on the whole phase plane of (1.3). 

For (5 = 0, we will have a system 

X = —y + {•y — a) X + bx^ — cx^ = P(x, y), 

(6.3) 

y = x = Q{x,y). 
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Applying the definition of a field rotation parameter to this system for the 
parameters a and c, respectively, we will get the following determinants: 



A, 



PQ'a - QP'a 



> 



(6.4) 



A, = PQ'^ - QP^ 



X* > 0. 



(6.5) 



It follows from (6.4) and (6.5) that on increasing a or c the vector field of (6.3) 
is rotated in positive direction (counterclockwise) on the whole phase plane 



For the study of multiple limit cycle bifurcations, we will use also two theorems 
by L.Perko (see [H]) which are formulated for the polynomial system (3.1). 

Theorem 6.1 (Wintner— Perko termination principle). Any one-para- 
meter family of multiplicity-m limit cycles of the relatively prime polynomial 
system (3.1) can he extended in a unique way to a maximal one-parameter 
family of multiplicity-m limit cycles of (3.1) which is either open or cyclic. 

If it is open, then it terminates either as the parameter or the limit cycles 
become unbounded; or, the family terminates either at a singular point of (3.1), 
which is typically a fine focus of multiplicity m, or on a (compound) separatrix 
cycle of (3.1), which is also typically of multiplicity m. 

The proof of this principle for the general polynomial system (3.1) with a vec- 
tor parameter /x G R" parallels the proof of the planar termination principle 
for the system 



with a scalar parameter A G R, since there is no loss of generality in assuming 
that system (3.1) is parameterized by a scalar parameter A (see |6], [H]). 

In particular, if A is a field rotation parameter of (6.6), the following Perko's 
theorem on monotonic families of multiple limit cycles is valid [IB] . 

Theorem 6.2. // Lq is a nonsingular multiple limit cycle of (6.6) for A = Aq, 
then Lo belongs to a one-parameter family of limit cycles of (6.6); furthermore: 

1) if the multiplicity of Lq is odd, then the family either expands or contracts 
monotonically as A increases through Aq; 

2) if the multiplicity of Lq is even, then Lq bifurcates into a stable and an 
unstable limit cycle as A varies from Aq in one sense and Lq disappears as A 
varies from Aq in the opposite sense; i. e., there is a fold bifurcation at Aq. 

Using these theorems and the results of the previous sections and applying 



of (6.3). 




(6.6) 
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the field rotation parameters of tlie systems (1-3) and (6.3), we will prove the 
following theorem. 

Theorem 6.3. The FitzHugh-Nagumo system (1.3) can have at most two 
limit cycles. 

Proof. First let us prove that system (1.3) can have at least two limit cycles 
supposing that all of the parameters of (1.3) are nonnegative (in the natural 
sense of the model). All other cases can be considered in a similar way. 

Let the parameters a, c, 7, 5 of (1.3) vanish and consider the quadratic system 

x = ~y + hx^, y = X. (6.7) 

This is a reversible system. It has a center at the origan O and cannot have 
limit cycles [6]. We will input the parameters 7, a, c, and 6 successively one 
by one into the system (6.7). 

Inputting a positive parameter 7, we will get a system 

X = —y + 'J X + bx'^, y = X, (6.8) 

the vector field of which is rotated in negative direction (clockwise) on the 
whole phase plane of (6.8). The origin of (6.8) becomes an unstable focus (or 
a node). 

Inputting a positive parameter a into (6.8), the vector field of the system 

X = —y + {'J — a) X + bx'^, y = x (6.9) 

will be rotated in positive direction (counterclockwise). For = 7, the origin 
becomes weak and changes the character of stability on further increasing a. 
The Andronov-Hopf bifurcation occurs for a = 7, and an unstable limit cycle, 
Fi, will appear from the origin [1], [6]. 

Inputting a positive parameter c into (6.9), we will get a cubic system 

X = —y + {'-y — a) X + bx'^ — cx^, y = x, (6.9) 

the vector field of which is also rotated in positive direction on the whole phase 
plane of (6.9). The structure and the character of stability of infinite singu- 
larities will be changed, and a stable limit, F2, surrounding Fi will appear 
immediately from infinity in this case. On further increasing the parameter c, 
the limit cycles Fi and F2 combine a semi-stable limit, F12, which then disap- 
pears in a "trajectory concentration" PP, 
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If to input a positive parameter S into (6.9), we will have again the original 
system (1.3). On further increasing 6, a saddle-node appears in the first quad- 
rant of the phase plane. It splits then in two singular points: a saddle S and 
an antisaddle A. Without loss of generality, we can fix the parameter 6, fixing 
the positions of the finite singularities O, S, A, and consider the system (1.3) 
with a positive parameter 7 which rotates the vector field of (1.3) on the whole 
phase plane. 

So, consider the original system (1.3) with a positive parameter 7. On increas- 
ing this parameter, the stable nodes O and A becomes first stable foci, then 
they change the character of their stability, becoming unstable foci. At these 
Andronov-Hopf bifurcations [1], [6], stable limit cycles will appear from the 
foci O and A. On further increasing 7, the limit cycles will expand and will 
disappear in small separatrix loops of the saddle S. If these loops are formed 
simultaneously, we will have a so-called eight-loop separatrix cycle. In this 
case, a big stable limit surrounding three singular points, O, S, and A, will 
appear from the eight-loop separatrix cycle after its destruction, expanding to 
infinity on increasing 7. If a small loop is formed earlier, for example, around 
the point O {A), then, on increasing 7, a big loop formed by two lower (upper) 
adjoining separatrices of the saddle 5* and surrounding the points O and A 
will appear. After its destruction, we will have simultaneously a big limit cycle 
surrounding three singular points, O, S, A, and a small limit cycle surrounding 
the point A (O). Thus, we have proved that system (1.3) can have at least 
two limit cycles, see also [i9\ for more detail. 

Let us prove now that this system has at most two limit cycles. The proof 
is carried out by contradiction applying Catastrophe Theory, see [6], [T8] . 
Consider the system (1.3) with three field rotation parameters: a, c, and 7 
(the parameters b and 6 can be fixed, since they do not generate limit cycles). 
Suppose that (1.3) has three limit cycles surrounding the only point, O, at the 
origin. Then we get into some domain of the parameters a c, and 7 being 
restricted by definite conditions on two other parameters, b and 6. This domain 
is bounded by two fold bifurcation surfaces forming a cusp bifurcation surface 
of multiplicity-three limit cycles in the space of the parameters a, c, and 7 

The corresponding maximal one-parameter family of multiplicity-three limit 
cycles cannot be cyclic, otherwise there will be at least one point corresponding 
to the limit cycle of multiplicity four (or even higher) in the parameter space. 
Extending the bifurcation curve of multiplicity-four limit cycles through this 
point and parameterizing the corresponding maximal one-parameter family 
of multiplicity-four limit cycles by a field rotation parameter, 7, according 
to Theorem 6.2, we will obtain two monotonic curves of, respectively, multi- 
plicity-three and one limit cycles which, by the Wintner-Perko termination 
principle (Theorem 6.1), terminate either at the point O or on an infinite 
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separatrix cycle surrounding this point. Since we know at least the cyclicity of 
the singular point which is equal to two (see [I9]), we have got a contradiction 
with the termination principle stating that the multiplicity of limit cycles 
cannot be higher than the multiplicity (cyclicity) of the singular point in 
which they terminate. 

If the maximal one-parameter family of multiplicity-four limit cycles is not 
cyclic, using the same principle (Theorem 6.1), this again contradicts the 
cyclicity of the origin (see [19]) not admitting the multiplicity of limit cy- 
cles to be higher than two. This contradiction completes the proof in the case 
of one singular point on the phase plane. 

Suppose that the system (1.3) has three finite singularities, O, S, A, and two 
small limit cycles around, e.g., the point O (the case when the limit cycles 
surround the point A is considered in a similar way). Then we get into some 
domain in the space of the parameters a, c, and 7 which is bounded by a fold 
bifurcation surface of multiplicity- two limit cycles [6], [T8] . 

The corresponding maximal one-parameter family of multiplicity-two limit 
cycles cannot be cyclic, otherwise there will be at least one point corresponding 
to the limit cycle of multiplicity three (or even higher) in the parameter space. 
Extending the bifurcation curve of multiplicity-three limit cycles through this 
point and parameterizing the corresponding maximal one-parameter family 
of multiplicity-three limit cycles by a field rotation parameter, 7, according 
to Theorem 6.2, we will obtain a monotonic curve which, by the Wintner- 
Perko termination principle (Theorem 6.1), terminates either at the point O 
or on some separatrix cycle surrounding this point. Since we know at least the 
cyclicity of the singular point which is equal to one in this case [TP], we have 
got a contradiction with the termination principle (Theorem 6.1). 

If the maximal one-parameter family of multiplicity-two limit cycles is not 
cyclic, using the same principle (Theorem 6.1), this again contradicts the 
cyclicity of O (see [TS]) not admitting the multiplicity of limit cycles higher 
than one. Moreover, it also follows from the termination principle that either 
an ordinary (small) separatrix loop or a big loop, or an eight-loop cannot have 
the multiplicity (cyclicity) higher than one in this case. Therefore, according 
to the same principle, there are no more than one limit cycle in the exterior 
domain surrounding all three finite singularities, O, S, and A. 

Thus, taking into account all other possibilities for limit cycle bifurcations (see 
[19j), we conclude that system (1.3) cannot have either a multiplicity-three 
limit cycle or more than two limit cycles in any configuration. The theorem is 
proved. □ 
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